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The Mott insulating phase of a one-dimensional bosonic gas trapped in optical lattices is described 
by a Bose-Hubbard model. A continuous unitary transformation is used to map this model onto 
an effective model conserving the number of elementary excitations. We obtain quantitative results 
for the kinetics and for the spectral weights of the low-energy excitations for a broad range of 
parameters in the insulating phase. By these results, recent Bragg spectroscopy experiments are 
explained. Evidence for a significant temperature of the order of the microscopic energy scales is 
found. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices provide tunable ex- 
perimental realizations of correlated many-particle sys- 
tems. Even the dimensionality of these systems can be 
chosen between one and three. Recently, the transition 
between superfluid and Mott insulating (MI) phases has 
been observed by tuning the ratio between the kinetic 
and the interaction energy 0, • Dynamical aspects can 
also be addressed experimentally 0, 13 • 

There are many theoretical treatments of these issues, 
ranging from mean-field approaches bp to more sophis- 
ticated techniques H |& 13 03 ET The transition 
at zero temperature (T = 0) has been analyzed in one- 
dimensional systems 0, H, 13 E3 Ell • But there is still a 
multitude of open questions. Our aim is to explain the 
dynamics of the excitations which has moved in the focus 
recently USUI. 

The present work is intended to clarify the significance 
of various excitation processes depending on the kinetic 
and the interaction energy and on the temperature. To 
this end, we will analyze the spectral weights of the one- 
dimensional MI phase with one boson per site n = 1. 

The low-energy excitations are either double occupan- 
cies (n — 2), which we will call 'particles', or they are 
holes (n = 0). Both excitations are gapped as long as 
the system is insulating. They become soft at the quan- 
tum critical point which separates the insulator from the 
superfluid. We use a continuous unitary transformation 
(CUT) m in to obtain an effective Hamiltonian in 
terms of the elementary excitations 'particle' and 'hole'. 
This effective Hamiltonian conserves the number of par- 
ticles and of holes 0] . The CUT is realized in real space 
in close analogy to the derivation of the generalized t- 
J model from the Hubbard model The strongly 

correlated many-boson problem is reduced to a problem 
involving a small number of particles and holes. This 



simplification enables us to calculate kinetic properties 
like the dispersions and spectral properties like spectral 
weights. 

The article is set-up as follows. First, the model and 
the relevant observable are introduced (Sect. II). Then, 
the method used is presented and described (Sect. III). 
The spectral weights are computed in Sect. IV. In Sect. 
V, the experiment is analyzed which requires calculations 
at finite temperatures as well. Finally, the results are 
discussed in Sect. VI. 



II. MODEL 

To be specific, we study the Bose-Hubbard model H = 
tH t + UHu in one dimension 

H = -tX>&+i + &t f i 6 i) + (U/2)Y,fM{n t - 1) (1) 

i i 

where the first term is the kinetic part tH t and the sec- 
ond term the repulsive interaction UHjj with U > 0. The 
bosonic annihilation (creation) operators are denoted by 

b^\ the number of bosons by hi = bib,. If needed the 
term = —fi hi is added to H to control the par- 
ticle number. For numerical simplicity, we truncate the 
local bosonic Hilbert space to four states. This does not 
change the relevant physics significantly 0, 13 @ • 

Besides the Hamiltonian we need to specify the excita- 
tion operator R. In the set-up of Refs.y and^the depth 
of the optical lattices is changed periodically to excite 
the system. In terms of the tight binding model JTJ this 
amounts to a periodic change of t and of U leading to 

R oc StH t + 5UH V . (2) 

Since multiples of H do not induce excitations we con- 
sider 

R^ R = R-{5U/U)H . (3) 
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Both operators R and R induce the same transitions 
(n\R\m) = (n\R\m) where |n) ^ |m) are eigen states 
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of H . Eventually, the relevant part of R is proportional 
to H t . For simplicity, we set the factor of proportionality 
to one. 

If the interaction dominates (U/t — > oo) the ground 
state of Q is the product state of precisely one boson 
per site |ref) = |l)i ® 1 1 ) 2 • •• <8> |1)jv, where \n)i denotes 
the local state at site i with n bosons. We take |ref) as our 
reference state; all deviations from |ref) are considered as 
elementary excitations. Since we restrict our calculation 
to four states per site we define three creation operators: 
h\\l)i — |0}i induces a hole at site i, p]\l)i — |2)j induces 
a particle at site i, and <4|l)i = |3)i induces a double- 
particle at site i. The operators h, p and d obey the 
commutation relations for hardcore bosons. 
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III. CUT 

The Hamiltonian Q conserves the number of bosons 
b. But if it is rewritten in terms of h,p, and d it is no 
longer particle-conserving, e.g., the application of H t to 
ref) generates a particle-hole pair. We use a CUT defined 
by 

d l H(l) = [r ) (l),H(l)} (4) 

to transform H(l = 0) = H from its form in Q to an 
effective Hamiltonian H c ff := H(l = 00) which conserves 
the number of elementary excitations, i.e., [Hjj, H e g] = 0. 
An appropriate choice of the infinitesimal generator r\ is 
defined by the matrix elements 

Vi,j(l) = sgn {q t - qj) H itj (l) (5) 

in an eigen basis of Hjj', q% is the corresponding eigen- 
value [14(. The structure of the resulting H c g and the 
implementation of the flow equation in second quantiza- 
tion in real space is described in detail in Refs. and 

The flow equation Q generates a proliferating num- 
ber of terms. For increasing I, more and more hardcore 
bosons are involved, e.g., annihilation and creation of 

three bosons plpl +1 p i+ 2PiPi + iPi + 2i an d processes over 
a larger and larger range occur, e.g., hopping over j 
sites p\ + jPi ■ To keep the number of terms finite we 
omit normal-ordered terms beyond a certain extension 
r. Normal-ordering means that the creation operators of 
the elementary excitations appear to the left of the anni- 
hilation operators. If terms appear which are not normal- 
ordered the commutation relations are applied to rewrite 
them in normal-ordered form. The normal-ordering is 
important since it ensures that only less important terms 
are omitted. Generically, such terms involve more parti- 
cles mull- 

We define the extension r as the distance between the 
rightmost and the leftmost creation or annihilation oper- 
ator in a term |l5l ITr| . The extension r of a term mea- 
sures the range of the physical process which is described 



FIG. 1: Residual off-diagonality for values oit/U £ {0.1,0.2} 
for maximal r — 2 and r = 3. Inset: Magnification for small 
values of the flow parameter I. The ROD for r = 3 and 
t/U = 0.2 displays non-monotonic behavior. In this case, we 
stop the flow at the first minimum of the ROD. Its position 
is indicated by a vertical line. 



by this term. So our restriction to a finite extension re- 
stricts the range of processes kept in the description. This 
is the most serious restriction in our approach. Note that 
the extension r implies for hardcore bosons that at max- 
imum r + 1 bosons are involved. 

For I < 00, H(l) still contains terms like p\h\ +1 which 
do not conserve the number of particles. To measure the 
extent to which H{1) deviates from the desired particle- 
conserving form, we introduce the residual off-diagonality 
(ROD) as the sum of the moduli squared of the coeffi- 
cients of all terms which change the number of elementary 
excitations. The ROD for maximal extensions r = 2 and 
r = 3 is shown in Fig. ^ Ideally, the ROD decreases 
monotonically with /. The calculation with r — 2 indeed 
shows the desired monotonic behavior of the ROD for all 
values oit/U. 

But non-monotonic behavior can also occur |16| . It is 
observed in Fig. for extension r = 3. While for small 
t/U the ROD still decays monotonically it displays non- 
monotonic behavior for larger t/U, see e.g. the ROD for 
t/U = 0.2 and r — 3 in Fig. 2] An uprise at larger 
values I }t 1/t signals that the intended transformation 
cannot be performed in a well-controlled way. The up- 
rise stems from matrix elements which do not decrease 
but increase, at least in an intermediate interval of the 
running variable I . Such an increase occurs if the number 
of excitations is not correlated with the total energy of 
the states. This means that two states are linked by an 
off-diagonal matrix element of which the state with low 
number of excitations is higher in total energy than the 
state with a higher number of excitations. The increase 
of such a matrix element bears the risk that the unavoid- 
able truncations imply a too severe approximation. 

The situation that the number of excitations is not cor- 
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related with the total energy can occur where the lower 
band edge of a continuum of more particles falls below the 
upper band edge of a continuum of less particles or be- 
low the dispersion of a single excitation, see, e.g., Ref.lT^. 
This situation implies additional life-time effects since the 
states with less particles can decay into the states with 
more particles. Even if the decay rates are very small the 
CUT can be spoilt by them because the CUT in its form 
defined by Eq. JSJ) correlates the particle number and the 
energy of the states. 

In order to proceed in a controlled way, we neglect the 
small life-time effects if they occur at all. This is done 
by stopping the CUT at Z m i n at the first minimum of the 
ROD. The position of l m i n for t/U = 0.2 and r = 3 is 
indicated in the inset of Fig. ^ by a vertical line. It is 
found that the remaining values of the ROD are small 
and thus negligible. Hence we omit the remaining off- 
diagonal terms. Only close to the critical point, where 
the MI phase vanishes, the present approach becomes 
insufficient. 
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FIG. 2: (color online) Phase diagram of the Mott insulating 
(MI) and the superfluid (SF) phase in the (t/U, fj,/U) plane. 
Dotted (dashed) lines show CUT results for maximal exten- 
sion r = 2 (r = 3). Solid lines (symbols) show series (DMRG) 
results from Refs. ITol and[3 



We check the reliability of our approach by comparing 
its results to those of other methods for the phase dia- 
gram in Fig. |21 The upper boundary of the MI phase is 
given by the particle gap A p := minfeW p (fc). The lower 
boundary is given by the negative hole gap — A h where 
A h := minfcu; h (fc). The dotted (dashed) curves result 
from CUTs for r = 2 (r = 3). For r = 2, the CUT can 
be performed till I — oo; for r = 3, the flow is stopped 
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FIG. 3: (color online) The dispersions of a particle ui p (k) 
(black) and a hole uj h (k) (green/grey) for t = 0.1U (solid), 
t = 0.1517 (dotted) and t = 0.2U (dashed). 



at / m in- Solid curves depict the findings by series expan- 
sion, the symbols those obtained by DMRG UGH- The 
agreement is very good in view of the truncation of the 
Hilbert space and in view of the low value of r. Note 
that the r = 3 result agrees better with the series and 
DMRG results than the r = 2 result. As expected, the 
deviations increase for larger values of t because longer- 
range processes become more important. Yet the values 
obtained by the CUT for the critical ratio x c := t/U, 
where the MI phase vanishes, are reasonable. We find 
xi r=2) = 0.271 and x { c r=3) = 0.258. By high accuracy 
density-matrix renormalization x c = 0.297 ± 0.01 was 
found, see Ref. [ll] and references therein. Series expan- 
sion provides x c = 0.26 ± 0.01 which is very close to our 

(r=3) 



value x c ~ ■ This fact underlines the similarity between 
series expansions and the real space CUT as employed in 
the present work. 

We conclude from the above findings for the phase dia- 
gram (Fig.[2J) that the mapping to the particle-conserving 
-ffeff works very well for large parts of the phase diagram. 
It does not, however, capture the Kosterlitz-Thouless na- 
ture of the transition itself [13. The CUT yields reliable 
results within the MI phase for t ^ 0.2U. Henceforth, re- 
sults for this regime will be shown which were obtained 
in the r = 3 scheme. 

In Fig. |31 the single-particle dispersion uj p (k) (one-hole 
dispersion w h (fc)) is shown as black (green or grey, resp.) 
curves. Both dispersions increase with t; the particle dis- 
persion always exceeds the hole dispersion. On increasing 
t, the center of the hole dispersion shifts to higher ener- 



gies while the center of the particle dispersion remains 
fairly constant in energy. 



IV. OBSERVABLE 

At the beginning of the flow (I — 0), the observable 
R is proportional to H t . The dynamic structure factor 
Sn(k, uj) encodes the response of the system to the appli- 
cation of the observable R. The observable R transfers no 
momentum because it is invariant with respect to transla- 
tions. Therefore, the Bragg spectroscopy 0,11] measures 
the response Sn(k = 0,u>) at momentum k = and en- 
ergy ui. A sketch of the spectral density Su(k — 0, w) 
for this observable is shown in Fig. We assume that 
the average energy is mainly determined by Hu in the 
MI regime. Then the first continuum is located at U 
and a second one at 2U. For small t/U the continua 
will be well separated. The energy-integrated spectral 
density in the first continuum is the spectral weight Si. 
Correspondingly, 5*2 stands for the spectral weight in the 
second continuum. 

To analyze the spectral weights in the effective model 
obtained by the CUT the observable R must be trans- 
formed as well. Before the CUT, the observable is 
R(l = 0) = Hf. It is a sum of local terms 

R(l = 0) = H t =J2 bLi /2 K + x/2 + bl+i/A-i/2 > (6) 
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FIG. 4: Sketch of the distribution of spectral weight S(k = 
0, uS). The weight centered around u) = U is given by Si, the 
weight around u) — 2U by 5*2- 
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where we have rewritten the sum as a sum over the bonds 
r. The bond positions are in the centers of two neighbor- 
ing sites. This notation emphasizes that the observable 
acts locally on bonds. The observable is transformed by 
the CUT to 



i? off = RQ = oo) = ^ R(l = oo, r) 



(7) 
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It is the sum over the transformed local observables R(l = 
oo, r) which are centered at the bonds r. 

The local observable is sketched schematically in Fig.[S] 
The sites on which the local observable acts are shown 
as filled circles. The state of the sites shown as empty 
circles is not altered by the observable. At I = 0, the 
observable is H t . It acts only locally on adjacent sites of 
the lattice as shown in Fig. [SJi. 

We compute the matrix elements of the excitation op- 
erator after the CUT, i.e., of R cS . This operator consists 
of terms of the form c (fi .„- l) (h^) n " {p r ) n *> ((fi) nd h mh p m "d md 
where we omitted all spatial subscripts denoting only the 
powers (n; m) = {nu n p ran m p ma) of the particular 
operator type. For the full expression we refer the reader 
to Appendix lAl The coefficient c^m) is the correspond- 
ing prefactor. The spectral weight 1^.^ is the integral 
of S(k = 0, uj) over all frequencies for momentum transfer 
k = 0. It stands for the excitation process starting from 
the states with rrid double-particles, m p particles, and rrih 



FIG. 5: Terms of the local observable R(l,r). The arrow 
indicates the bond on which a term acts. Filled circles: sites 
on which a non-trivial operator acts; 'non-trivial' means that 
the operator is different from the identity. Open circles: sites 
on which the term under study does not act, i.e., the term acts 
as identity on these sites, (a) At I = the observable is Ht- 
It is composed of local terms that act only on adjacent sites 
of the lattice, (b) and (c) More complicated terms appear 
during the flow. The sum of the distances of all operators in 
a term is our measure ro for the extension of a term. It is 2 
for (b) and 4 for (c). Terms beyond a certain extension are 
omitted. 



holes and leading to the states with rid double-particles, 
n p particles, and holes. 

In the course of the flow, contributions to the observ- 
able appear which do not act on the bond on which the 
initial observable is centered. The initial local process 
spreads out in real space for I — > oo. Examples are 
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sketched in Fig. [SJa-c. In order to avoid proliferation, 
also the terms in the observable have to be truncated. 
Like for the terms in the Hamiltonian we introduce a 
measure for the extension of the terms in the observable. 
This measure, however, is slightly different: It is the sum 
ro of the distances of all its local creation or annihilation 
operators to r. If the value tq of a certain term exceeds 
a preset truncation criterion, this term is neglected. 

Figures ]Ejp-c illustrate the truncation criterion for the 
observable. An operator that acts on the sites shown in 
Figure 03 meets the truncation criterion for tq = 3. It is 
kept in a rg = 3 calculation. An operator acting on the 
sites in Figure is kept in the calculation with ro = 4. 
But it does not meet the truncation criterion for tq = 3; 
in a calculation with ro = 3 it is discarded. 

Once the local effective observables R(l = oo, r) are 
known, the total effective expression R cS is given by the 
sum over all bonds Q. 

The results for the observable truncations ro £ 
{2,3,4} are shown in Fig.H3 At zero temperature, only 
the processes starting from the ground state are relevant. 
These processes are the (n; 0) processes which start from 
the reference state |ref) because the CUT is constructed 
such that the reference state, which is the vacuum of ex- 
citations, becomes the ground state after the transforma- 
tion [lj, IS ■ The weights of the (n;0) processes are 
shown in Fig. The by far dominant weight is ^( C jf p . ) ; 

the process ^y^p-o) * s l° wer by orders of magnitude. The 
agreement between results for various truncations is very 
good. 

The particle-hole pair excited in If^ p .m has about the 
energy U for low values of t/U, while If^p-o) ^ ea -ds to a 
response at 2U. Hence we find noticeable weight only 
around U in accordance with recent quantum Monte 
Carlo (QMC) data [13 . 

At finite temperature, excitations are present before 
i? off is applied. At not too high temperatures, inde- 
pendent particles p^ and holes h' are the prevailing ex- 
citations. Other excitations, for instance d< or corre- 
lated states (p^) 2 of two p-particles, are higher in en- 
ergy and thus much less likely. So the processes starting 
from a particle or a hole are the important ones which 
come into play for T > 0. Hence we focus on i^ p . p j, 

^(°2h P -h)' an d I(hd-p)- These weights are shown in Fig. EJ>. 
The results for ij^. p j depend only very little on trunca- 
tion. A larger dependence on the truncation is found for 
■^(h2p-p) +^°2h P -h)- the agreement is still satisfactory. 

The processes I^ 2p -p) an< ^ ^(2hp-h) increase the energy 
by about U because they create an additional particle- 
hole pair. The InJd-p) P rocess increments the energy by 
about 2U because a hole and a double-particle is gener- 
ated. 
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FIG. 6: (color online) Spectral weights of various processes 
for truncations ro £ {2,3,4}. (a) I?h p; o) (black curves) and 
^( e 2 ff h2p;0) (green/grey), (b) J ( f 2p;p) +/ ( f hp;h) (black) and J ( f d;p) 
(green/grey) 



V. APPROACHING EXPERIMENT 



Let us address the question what causes the high en- 
ergy peak in Refs. 0,0,9- It was suggested that certain 
defects, namely an adjacent pair of a singly and a doubly 
occupied site, are at the origin of the high energy peak. 
The weight of such processes for a given particle state is 
quantified by ^(hd P ) ■ ^s relatively high value, see Fig.Et, 
puts the presumption that such defects cause the peak at 
2U on a quantitative basis. 
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A. Zero Temperature 

But what generates such defects? At zero or at very 
low temperature, the inhomogeneity of the parabolic trap 
can imply the existence of plateaus of various occu pat ions 
(n) G {0, 1,2,.. .} depending on the total filling [TJ, US 
|20|. In the transition region from one integer value of 
(n) to the next, defects occur which lead to excitations 
at 2U. 

Yet it is unlikely that this mechanism explains the ex- 
perimental finding since the transition region is fairly 
short at low value of T/U and t/U, i.e., the plateaus 
prevail. The high energy peak at 2U, however, has less 
weight by only a factor 2 to 5 compared to the weight in 
the low energy peak at U |3j. So we conclude that the 
inhomogeneity of the traps alone cannot be sufficient to 
account for the experimental findings. 




T/U 



FIG. 7: (color online) Temperature dependence of the relative 
spectral weight S2/S1 for various values of t/U. Inset: S2/S1 
as function of t/U for various temperatures. 



B. Finite Temperature 

At higher temperatures, thermal fluctuations are a 
likely candidate for the origin of the defects. Hence we 
estimate their effect in the following way. Thermally in- 
duced triply occupied d states are neglected because they 
are very rare due to their high energy of ui d (k) — 2/x » 2U 
above the vacuum |ref). We focus on the particle p and 
the holes h. The average occupation of these states is 
estimated by a previously introduced approximate hard- 
core boson statistics [2l| 

«) = exp {-P{w'{k) - n) (n vac ) (8a) 
(n vac ) = (l + zP^ + A/?))^ , (8b) 

where z a = (27T)- 1 dke~^ u<r ^-^\ fx = /i p = -^ h . 
The chemical potential changes sign for the holes be- 
cause a hole stands for the absence of one of the original 
bosons b. Equation JSJ is obtained from the statistics 
of bosons without any interaction by correcting globally 
for the overcounting of states. The fact that the over- 
counting is remedied on an average level implies that® 
represents only an approximate, classical description |22| . 
The chemical potential fi is determined self-consistently 
such that as many particles as holes are excited, i.e., the 
average number of bosons b per site remains one. 

Next, we turn to the computation of spectral weights 
at finite temperatures. At not too high temperatures, 
the relevant channels are (hp;0), (h2p;p), (2hp;h), which 
excite at around U and (2h2p;0), (hd;p), which excite 
at around 2U. The spectral weights at finite T are cal- 
culated by Fermis golden rule as before. The essential 
additional ingredient is the probability to find a parti- 
cle, or a hole, respectively, to start from and to find sites 
which can be excited, i.e., which are not occupied by a 



particle or a hole. This leads us to the equations 

Toff,T _ Tcft /„vac\n h +n p +n d /q \ 

1 (n;0) ~ J (n;0)\ n / l ya ) 

/ vac\n h +np+n d r2ir 

C) = JLA ^ — I #?> ( 9b ) 

with a G {p,h}. The powers of (n vac ) account for the 
probability that the sites to be excited are not blocked by 
other excitations; the factor (ro£) accounts for the prob- 
ability that the excitation necessary for the particular 
process is present. The momentum dependence in I k < cS 
stems from the momentum dependence of the annihilated 
particle or hole. It is computed by the sum of the moduli 
squared of the Fourier transform of the matrix elements 
C(n;m) in the real space coordinate of the annihilated ex- 
citation. 

The spectral weight in the high energy peak (at 2U) 
relative to the weight in the low energy peak (at U) is 
given by the ratio S 2 /Si where S 2 = ^ 2 h2 Pi o) + 7 (hd; P ) 
and Si = /(fpfo) + /gp. p) + / ( ° 2 ff hp;h)- Fig. □ displays this 
key quantity as function oit/U and T/U for tq = 4. The 
difference to the result for tq — 3 is less than 0.005. 

Two regimes can be distinguished. For low values of 
T ^ 0.19U, the ratio S2/S1 increases on increasing t/U 
because the particle gap A p decreases so that (n-^~ ) 
grows. For higher values of T ^ 0.19C7, the increase 
of («^~o) is overcompensated by the decrease of the 
weights ^(2h2p-o) + ^(hd^>) an d t ne increase of the weights 

J S) + 7 Sp) + 7 Sp;h)' cf - Fi S- El so that the rela- 
tive spectral weight decreases on increasing t/U . Around 
T = 0.19C/, the ratio is fairly independent of t/U . 

The experimental value of S2/S1 J| is about 0.2 — 
0.5 for small values of t/U (t < Q.03U). It increases on 
approaching the superfluid phase. Hence, our estimate 
for S2/S1 implies a significant temperature T w U/3 in 
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the MI phase, which was not expected. This is the main 
result of our analysis of the spectral weights. 



VI. DISCUSSION 

The analysis of the spectral weights lead us to the 
conclusion that the temperature of the Mott-insulating 
phases must be quite considerable. At first sight, this 
comes as a surprise because other experiments on cold 
atoms imply very low temperatures, see for instance Ref. 
2. But the seeming contradiction can be resolved. 

One has to consider the entropies involved |U |24| . 
The entropy per boson S/N in a three-dimensional har- 
monic trap is given by « 3.6(1 — /o) where /o is the 
fraction of the Bose-Einstein condensate j^. Adia- 
batic loading into the optical lattice keeps S/N con- 
stant so that we can estimate the temperature of the 
MI phase from the derivative of the free energy [2l| . 
For fo = 0.95; 0.9; 0.8 we obtain significant tempera- 
tures T/U = 0.12;0.17;0.27in satisfactory agreement 
with the analysis of the spectral weights [22| • The anal- 
ogous estimate in the case of the Tonks-Girardeau limit 
(large U and average filling n below unity, see Ref. 0) 
leads to temperatures of the order of the hopping t: 
T/t = 0.17; 0.32; 0.61 (assuming n = 1/2) H^. These 
values are again in good agreement with experiment. 

The physical interpretation is the following. On ap- 
proaching the Mott insulator by changing the filling to 
the commensurate value n — > 1 the temperature has to 
rise because the available phase space decreases. For 
n < 1 there is phase space without occupying any site 
with two or more bosons b because one can choose be- 
tween occupation or 1. But at n = 1 the state without 
any doubly occupied site is unique. Hence no entropy 
can be generated without inducing doubly occupied and 
empty sites, i.e., excitations of p- and of h-type. This in 
turn requires that the temperature is of the order of the 
gap which agrees well with our analysis of the spectral 
weights. 

We like to draw the reader's attention to the fact that 
our result of a fairly large temperature (T ~ U/3) pro- 
vides also a possible explanation why no response at 
« 2U was found by QMC a low temperatures [12|] Fur- 
ther investigations will be certainly fruitful, e.g., it would 
be interesting to obtain QMC results for the excitation 
operator R that we used in o ur p resent work. The ex- 
citation operator used in Ref. [12 vanishes for vanishing 
momentum so that it is less suited to describe the exper- 
imental Bragg spectroscopy |3J. 

In the attempt to provide quantitive numbers for the 
temperature in the MI bose systems one must be cautious 
because the bosonic systems are shaken fairly strongly in 
experiment H0. Though it was ascertained that the 
experiment was conducted in the linear regime it might 
be that the systems are heated by the probe procedure so 
that the temperatures seen in the spectroscopic investi- 
gations are higher than those seen by other experimental 



investigations. 

It would be rewarding to clarify the precision of the 
spectroscopic investigations. Then the pronounced T de- 
pendence of the relative spectral weight in Fig. [7| could 
be used as a thermometer for Mott-insulating bosons in 
optical lattices. 

In summary, we have studied spectral properties of 
bosons in one-dimensional optical lattices using particle- 
conserving continuous unitary transformations (CUTs). 
At T = and for small t/U spectral weight is only present 
at energies ks U . Recent experimental peaks at f=a 2Z7 Q 
can be explained assuming T sa U/3. Our results suggest 
to investigate the effects of finite T on bosons in optical 
lattices much more thoroughly. 
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APPENDIX A: EXPLICIT EXPRESSION FOR 
THE OBSERVABLE 

The local observable depending on the flow parameter 
I reads 



R{l,v) 



Y, c(i,{i,i',j,f,k,k'}) 

{i,i',j,j',k,k'} 



ii+r 



h L+rh+r 



h i' m +r X 



P ]i+r ' ■■• ' Pj np +rPj[+r ' — ' Pj'+r X 



' d kn d +r d k' 1 +r 



■d h , 



k' +r- 



(Al) 



The numbers (rih n p nd',rn,h m p md) are defined as the 
number of operators involved in a term in Eq. IA1I The 
number of creation operators h\ p\ and dJ is given by 
n/j, n p , and n^, respectively. The number of annihilation 
operators h, p, and d is given by m/,, m p , and mj, respec- 
tively. A set of these six numbers defines the type (n; fh) 
of a process. The variables {i, i', j, j', k, k'} are multi- 



indices, e. g. i = {ii, 



,} which give the position of 



the operator. The coefficients c(l,{i,i',j,j',k,k'}) keep 
track of the amplitudes of these processes during the flow. 
Their value at I = oo defines the effective observable R eS . 
The total observable is the sum 



i?(0 = ]>>(?, r). 



(A2) 



No phase factors occur so that the observable is invariant 
with respect to translations. Hence no momentum trans- 
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fer takes place and the response at k — is the relevant one. 
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